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The dynamics of a two-component dilute Bose gas of atoms 
at zero temperature is described in the mean field approxima- 
tion by a two-component Gross-Pitaevskii Equation. We solve 
this equation assuming a Gaussian shape for the wavefunc- 
tion, where the free parameters of the trial wavefunction are 
determined using a moment method. We derive equilibrium 
states and the phase diagrams for the stability for positive 
and negative s-wave scattering lengths, and obtain the low 
energy excitation frequencies corresponding to the collective 
motion of the two Bose condensates. 



The paper is organized as follows. In Sect. 2 of this 
paper we first define our model, and derive the equations 
of motion for the parameters of our Gaussian wave func- 
tion. In particular, we study in detail the simple case of 
of isotropic traps and isotropic condensates. In Sec. 3 
we investigate first the case where the centers of the two 
trapping potentials for the two condensates coincide. For 
this configuration we calculate the equilibrium points and 
the low energy eigenmodes. In addition, we analyze the 
stability of the system. Finally, in Sec. 4 we consider 
the case of displaced trap centers, and compare with the 
results of Sec. 3. 
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I. INTRODUCTION 

Since the first observations of Bose-Einstein-conden- 
sation in an alkali vapor considerable effort has 

been made to characterize these systems both experi- 
mentally and theoretically. One of the latest advances 
is the experimental preparation of two-component Bose- 
condensates in different internal states, exhibiting a spa- 
tial overlap Q| . In this experiment ^''Rb atoms in two dif- 
ferent hyperfine states, |1) — |1, —1) and |2) = |2,2), are 
loaded into a magneto-optical trap, are cooled down to 
effectively zero temperature and subsequently undergo a 
phase transition towards a Bose-Einstein-condensate . 

The Gross-Pitaevskii Equation (GPE), a nonlinear 
Schrodinger Equation (NLSE) for the macroscopic wave 
function of the Bose condensed gas, provides an accu- 
rate description of the ground state and of the excitation 
spectrum of a dilute Bose condensate at zero temperature 
|p|-pO[ . Recently, several groups have solved the GPE for 
single component Bose condensates and found excellent 
agreement with experiment jl^ . The aim of this pa- 
per is to study solutions of the two- component GPE for a 
harmonic trapping potential. In particular we will derive 
the Bose ground state, investigate its stability proper- 
ties, and determine the low-lying excitation frequencies 
for positive and negative scattering lengths. Our tech- 
nique of solving the two-component GPE is based on a 
Gaussian ansatz for the condensate wave function where 
the open parameters are determined by moment methods 
pl[ (which is equivalent to a variational technique [Q). 
This allows one to obtain essentially analytical solutions 
of the GPE, and thus complements numerical studies of 
this equation pJ. 



II. BASIC EQUATIONS 

A. The model 

Zero-temperature mean-field theory provides an accu- 
rate theoretical description of dilute Bose condensed sys- 
tems. In particular, the dynamics of a two-component 
Bose gas can be modeled by the coupled Gross-Pitaeves- 
kii equations 
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which are a set of nonlinear Schrodinger equations for the 
macroscopic wave functions 5'i and ^'2 for the two com- 
ponents of the condensate. In Eqs. (|l|) Vn,{n = 1,2) de- 
notes the trap potentials which we assume to be harmonic 
in agreement with common experimental situations 
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K(r) = -Mlo^ {\l{x - xn? + Xy + \l7?) 



(2) 



The parameters Ar; (j] — x, y, z) account for the anisotro- 
py of the trap, and the trap centers for the first and 
second component are displaced in the x-direction by 
xi and X2, respectively. The coupling constants Unm 
are related to the scattering lengths a„„i by Unm = 
inh^anm/M. In this paper we allow only for the follow- 
ing elastic scattering processes: |1)|1) llljl)' |2)|2) — > 
1 2) 1 2), |1)|2) 1 1)1 2). Up to now, Eqs. (|j) have been 
solved in the Thomas-Fermi approximation, exploring 
the spatial densities of the mixtures |22] and the long 
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wavelength excitations [g^. Numerical studies of these 
systems have been carried out in [^^ . 

For non-interacting particles the Eqs. (|l|) are solved 
for the ground state by Gaussian functions, which are no 
more solutions in the general case due to the nonlinear 
self-interaction and coupling terms. The basic assump- 
tion behind the following derivations is that, in the case of 
interactions, the wave function is still well approximated 
by a Gaussian, 

*„(f,i) = A„(i)e-M/3-(*)(--"o)^+/3„„(t)y^+/3„.(0-^] (3) 

where the (3nr] — Pnri + complex numbers. The 

adjustable parameters in it can be interpreted as the am- 
plitudes A„, the width w„^ = 1/ ^ Pn-q ^'^'^ the curvature 

(M„^y^)-i/2 = {Pl^y^/^ of the wavefunction. To 

derive reliable results for the corresponding equations of 
motion it is necessary to take also the imaginary part 
of the wavefunction into account as shown in j25| for a 
similar case. One has to keep in mind when applying 
this ansatz to problems where losses are included, the 
loss terms usually induce aberrations, thus making the 
profile non-Gaussian. When nonlinear loss or gain terms 
are substantially important these aberrations make the 
information that can be obtained form an ansatz like Eqs. 

only qualitative. However, in the present case when 
losses are small corrections and do not determine the dy- 
namics the solutions of Eq. can be well approximated 
by Gaussian functions. 

In the following we find it convenient to work with 
scaled variables. We will measure energies in units of 
fiuj^ and will scale lengths with respect to ao\/2 where 
ao = {Mu)]^/"^ is the size of the ground state of the 
bare harmonic oscillator. 



B. Moment equations 

In Ref. solution of the single component Gross- 
Pitaevskii equation was formulated as a variational prob- 
lem, where a Lagrangian density C is derived so that 
the stationary point of the corresponding action gives 
the nonlinear Schrodinger equation. An approximate so- 
lution of the GPE was then derived by finding the ex- 
tremum of the action within a set of Gaussian trial func- 
tions which give "Newton type" equations for the varia- 
tional parameters. In contrast to such a variational ana- 
lysis , we will solve the two-component GPE by means 
of a moment-method. In the case of real Unm, i.e. no 
losses this method is equivalent to the variational 
approach. In principle, this approach is more general, as 
it can also be applied in a situation, where loss terms are 
included in Eqs. (Q). 

Using the Gaussian ansatz for a system where one con- 
densate is centered at xi = and the second one at 
X2 = a it is possible to calculate the number of particles 



in each condensate n = 1 , 2 as a function of the parame- 
ters 



dV|*„(f,i)p = |^„(i)p7r3/2w„^w„j,w 



(4) 



In a similar way, the second moments of the spatial co- 
ordinates rj ^ x,y, z and the second moments of the mo- 
menta are given by 
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respectively. These fourteen variables provide a complete 
description of the system within the previously men- 
tioned approximation. To obtain a closed set of equations 
of motion for A'„(t), w„„(i) and Mnnit), we take the time 
derivative of Eqs. (§,||,^, eliminate the time derivatives 
of the wave function with the help of the Schrodinger 
equations (|l|), and work out the resulting integrals over 
r with the help of the Gaussian ansatz (m . We obtain 
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where we replaced qnm = UnmVSn^ The corre- 

sponding equations for the second condensate are ob- 
tained by exchanging every index 1 and 2 and the cor- 
responding equations for the z-direction are obtained by 
exchanging the indices y and z in the equations for the 
y-direction. 



III. ISOTROPIC TRAPS AND IDENTICAL 
CONDENSATES: NO DISPLACEMENT OF TRAP 
CENTERS 



A. Equations 



trapping potential, the second term results from the ki- 
netic energy term, and the third and fourth term are due 
to the interaction between alike particles and between 
particles in different states. Below we further reduce the 
number of open parameters by assuming the scattering 
lengths ail and 022 to be equal. 



B. Equilibrium Points 

If we now examine the state in equilibrium of the sys- 
tem, we find that the requirement of equal scattering 
lengths implies that the condensates must have equal 
widths wi = W2 = w. Equilibrium points are obtained 
from the extrema of the potential (||), which gives 
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In the hydrodynamic approximation this equation of 
motion reduces further and can be analytically solved 
as w^g" = A^^/^(gii + 912)^/^- To test for stabili- 
ty of these equilibrium points requires us to check 
for minima with the help of the Hessian-matrix H — 
[d^Ves/ {dwindw2n)]w^q- This analysis will be carried out 
in the following subsection. 



The simplest case which is amenable to analytical 
treatment corresponds to an isotropic trapping potential 
Arc = Aj, = Az = 1 with (non-displaced) trap centers Xn — 
(n = 1,2) We will study the case of a fixed number 
of particles in each condensate, U(i = U22 = UI2 = 0, 
and, in particular, consider them equal Ni = N2 — N . 
With these assumptions, and the initial condition of two 
isotropic condensates Wnx = w„j, = Wnz = w„ (n = 1,2), 
the equations of motion for the widths of the condensates 
have the form 
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whereby exchanging the indices we get the corresponding 
one for W2. 

Similar to Ref. we can find a "potential" for these 
equations. Hence, the evolution of the width of the con- 
densates can be viewed as the coordinates of a fictitious 
classical particle moving in a three dimensional potential 
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C. Low energy eigenmodes and Stability 

To obtain the frequencies of small amplitude oscilla- 
tions around the equilibrium positions we linearize the 



full equations (H) around the equilibrium points. 



Swnrj- We find the following four frequencies: 
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which can be interpreted as collective oscillations de- 
picted in Fig. 1. The modes (a) and (6) are two fold 
degenerate, corresponding to oscillation e.g. in the x — y- 
plane and the x — z-plane 1 28 1 . Modes (c) and [d) corres- 
pond to isotropic oscillations. 



where the first term on the RHS stems from the harmonic 
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FIG. 1. Diagrammatical representation in the x-y-plane of 
the four collective excitations with the lowest energy. Oscilla- 
tions of particles in state |1) are visualized by solid line arrows, 
oscillations of particles in state |2) by dashed line arrows. 



The order in which these motions appear with increas- 
ing energy depends strongly on the values of the scatter- 
ing lengths involved. In Fig. 2 we display the spectrum 
for a positive value of an. For large, negative values of 
ai2, the two modes with the lowest energies are those 
with the maximum spatial overlap, namely mode b and 
mode c. If we now consider the region where ai2 is posi- 
tive and large, we naturally find that the modes a and d, 
for which the spatial overlap is minimized, are the lowest 
energy modes. The point at which this crossing happens 
changes with changing values of an. Again we find, as 
in [Q, that the widths of the condensates remain finite 
up to the point of collapse. 
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FIG. 2. Spectrum of the low lying excitation frequencies for 
a fixed, positive value qn — 2.3 x 10^**. With qi2 increasing 
from negative values to positive ones the two low lying modes 
and the two high lying ones get exchanged. 



The eigenfrequencies (O) provide the stability region 



of the equilibrium pointspo]) 
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Since for ai2 = we recover the single component results 
[p^ where only particles in one state are present, we see 
that for positive values of ai2 the region of stability in- 
creases and for negative ones it decreases. This is due to 
contributions to the total system energy from either the 
repulsive or attractive interaction. However, this effect 
is non-linear due to the nature of this interaction. It is 
illustrated in Fig. 3. 
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FIG. 3. Phase diagram for the stability of the two-condensate 
system for — 2000 particles. For positive ai2 the region of 
stability is increased due to the positive energy contribution from 
the repulsive interaction. However, it does not expand linearly 
with increasing ai2. 
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IV. ISOTROPIC TRAPS AND IDENTICAL 
CONDENSATES: DISPLACED TRAP CENTERS 

A. Equations 

Let us now consider the case of two spatially displaced 
traps and work out the differences to the above situation. 
We assume a displacement in the x-direction and make 
an ansatz analogous to Eqs.(||) and (||) where one con- 
densate is centered at xi = and the second at X2 = a. 
Also in this case we assume that the number of parti- 
cles is conserved, what implies that the imaginary parts 
of the coupling constants are again zero. We apply the 
moment method as in the previous case and arrive, since 
the spatial symmetry is now broken in the a;-direction, at 
the following set of equations of motion for the widths of 
the two condensates in the x-direction, w„a,, and in the 
y- and z-directions, w„ 

1 iViqii 
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Again one gets the equations for W2x and W2 by simply 
exchanging the indices. We can also like, in the non- 
displaced case, find a potential of the form of Eq.(||), 
which is suitable for the derivation of these equations. 
However, in this case the term originating from the inter- 
action of the particles in the different states is multiplied 



by an additional factor, exp[— q;^/(wJ^ 



''2x 



)] , which en- 



sures that the influence of these interactions gets weaker, 
the more the condensates are displaced. This then im- 
plies an even smaller region of spatial overlap. Note that 
the widths of the condensates are included in this factor 
and it therefore has a non-negligible effect on the physics 
of the system in addition to weakening the magnitude of 
effects found in the non-displaced case. 
Let us start calculating the equilibrium states for this 
case from Eqs.(|l|). In Fig. 4 we plotted the widths Wx 
and w of the condensates for an and ai2 positive and 
equal. Once the condensates are separated by a distance 
a — \J (w^j, -|- W22;)/2 the repulsive interaction between 
particles from different states becomes an effectively at- 
tractive one (and v.v.), as can be seen from the fact that 
the width on the x-direction gets smaller then the width 
in the y- and z-direction. The explanation for this effect 
is, that the direction of the repulsive force with respect to 
the surfaces of the condensates changes in these regions 
where they do no longer overlap. At large distances when 
the interaction between the particles in different states 
becomes negligible again, the case of two independent, 
isotropic condensates is recovered again. 
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1 2 3 a 4 5 6 

FIG. 4. Width of the condensates in the a;-direction and in 
the y-,2:-directions for qh = 0.2 x 10* and qvi — 0.2 x 10*. 
Due to the repulsive interaction between particles from differ- 
ent states, the condensate gets squeezed in the a;-direction at 
a certain distance and then expands once the region of overlap 
diminishes further. 



B. Low energy eigenmodes and Stability 

Looking at the excitation spectrum of the two compo- 
nent system by analyzing it again via a linear expansion 
around the equilibrium points, we find, as can be seen 
from Fig. 5, that the degeneracies existing in the non- 
displaced case get lifted with increasing distance. Bring- 
ing the condensates that far from each other that they 
can be considered as two independent condensates |lj] 
the spectrum reduces again to two three-fold degenerate 
values, since the condensates become again isotropic. 




1 2 3 (x 4 5 6 

FIG. 5. Spectrum of the low lying excitation energies for 
Qii = q\2 = 6 X 10^*. The degeneracies are lifted with in- 
creasing a and for o > 2 it reduces to the case of two .isotropic, 
single condenstates again. 



In this Gaussian-ansatz we can again easily interpret 
these modes to be of the following types: 
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(14) 



The members of these vectors represent the oscillation- 
amphtude in the directions {5ix 5iy 5iz ^2x ^2y ^2z), 
where 7; (/ = a, 6, c, d) depends on the parameters of the 
system. The hfting of the degeneracies is mainly due to 
the fact, that the spatial displaced system is not invariant 
under an exchange of the variables wi^ and wi^ and not 
of vf2y and W22, as can be seen from Eqs. ([l3|). 

Let us now look at the stability of the stationary states, 
which we again retrieve from the analysis of the Hes- 
sian matrix. As can be seen from Fig. 6 the further the 
condensates separate from each other, the weaker the ef- 
fect of additional stabilization (destabilization) becomes, 
since it results from the repulsive (attractive) interaction 
between the condensates. For a 1.5 one recovers the 
stability condition for the case of two separate conden- 
sates. 
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FIG. 6. Stability region for displaced traps with a displace- 
ment parameter a. For far displaced traps the result for inde- 
pendent condensates is recovered. 



V. CONCLUSIONS 

We used a moment method to describe the behavior 
of two overlapping Bose-Einstein-Condensates in differ- 
ent internal states, assuming that the shapes of the in- 
dividual wave functions do not deviate too much from a 
Gaussian-shape. With this method it is straightforward 
to derive equations of motion for the open parameters 
defining the exact form of the Gaussian function, espe- 
cially the width of the condensates. They can for the 
general cases easily treated numerically and for special 
cases even analytical solutions are possible. 

For the case of isotropic traps and isotropic conden- 
sates we did calculate the equilibrium points for the sys- 



tem and analyzed their stability. We found, that in com- 
parison to the case of one pure condensate the region 
of stability can be extended (shrinked) if the interac- 
tion energy between particles in different states is pos- 
itive (negative). Due to the nonlinear character of the 
interaction also this effect is of nonlinear nature. It ap- 
pears strongest for non-displaced traps and gets weaker 
the more the centers of the two traps are apart from each 
other. We were also able to calculate the energies of the 
first excited states and interpret the collective motions. 
For the case of displaced traps the degeneracies existing 
in the non-displaced case get lifted. 
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